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METHOD AND APPARATUS FOR ALLOCATING LINK 
BANDWIDTH AS FUNCTION OF QOS REQUIREMENT 

RELATED APPLICATION 

This application claims priority pursuant to 35 U.S.C. Sec 1 19(e) to US 
Provisional Application No.60/455,030, filed March 13, 2003, entitled METHOD 
AND APPARATUS FOR ALLOCATING BANDWIDTH TO A LINK OF A 
5 COMMUNICATIONS NETWORK AS FUNCTION OF QOS REQUIREMENT, 
the subject matter thereof being fully incorporated by reference herein. 

FIELD OF THE INVENTION 

The invention relates generally to network resource allocation methods, and 
particularly to the allocation of bandwidth for links of a communications network. 

10 BACKGROUND OF THE INVENTION 

A fundamental problem of Internet traffic engineering is that of determining 
the bandwidth (bits/sec), or link speed, required to carry a traffic load (bits/sec) 
offered to a given transmission link, and to satisfy specified quality-of-service 
(QoS) requirements for the offered traffic. The offered load is comprised of packets 

15 of varying sizes arriving for transmission on the link. Arriving packets can queue 
up at an input buffer for the link, and are dropped if the queue size (in bits) is larger 
than the input buffer size (in bits). Exemplary QoS metrics for the arriving packets 
are queuing delay and packet loss. 

The bandwidth required to satisfy a specified QoS requirement (hereafter 

20 sometimes referred to as the "QoS bandwidth") depends on the queue-length 

process which, in turn, depends heavily on the statistical properties of the packet 
arrivals and sizes. These statistical properties also change as the mean connection 
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load changes. Because of this statistical variation, the bandwidth allocation 
methods of the art have often produced sub-optimal bandwidth allocation. 

SUMMARY OF THE INVENTION 

A method is provided for determining a bandwidth allocation to provide a 
5 specified QOS requirement that takes appropriate account of statistical variations for 
packet streams in a transmission link. In particular, the method of the invention 
forms a statistical model of the packet stream using fractional sum difference 
statistical models and evaluates that model in respect to synthetically generated 
traffic streams. The bandwidth allocation approach of the invention is specified in 
10 terms of the bandwidth, p, required for a traffic load, x, subject to the requirements 
of a maximum queuing delay, 5, and a packet loss limitation parameter, co . 
Accordingly, that bandwidth allocation approach is implemented as a statistical 
model for p as a function of x, 5 and co . 

In a further embodiment of the invention, an algorithm is determined from 
15 the statistical model so developed for specifying the QOS bandwidth needed in a 
given application and for a given transmission link. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 provides a flow chart directed to the bandwidth allocation method 
of the invention. 

20 Figure 2 shows a data distribution for modeling a bandwidth allocation 

model according to the invention. 

Figure 3 shows a data distribution for modeling a refinement in that 
bandwidth allocation model. 
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Figure 4 shows a data distribution in respect to a further refinement in that 
bandwidth allocation model. 

Figure 5 shows a data distribution directed to a specification of a statistical 
parameter term for a bandwidth allocation model of the invention. 
5 Figure 6 shows a plot of QoS utilization for a bandwidth allocation model of 

the invention. 

DETAILED DESCRIPTION OF THE INVENTION 

The invention disclosed herein is directed to bandwidth allocation for links 
of a communication network. In particular, the invention provides an improved 

10 methodology for determining what link bandwidth is needed for a given traffic bit 
rate to ensure that the transmission on the link maintains a defined QoS. It is noted 
that the inventors have authored a companion article to the disclosure of the 
invention herein, which may be useful as a source of background detail or for 
additional information respecting the fractional sum difference statistical models 

15 used in the development of the invention (as discussed more fully below). That 
article, entitled "A Statistical Model For Allocating Bandwidth To Best-Effort 
Internet Traffic," 2004 [In Preparation] is hereby incorporated by reference herein. 

Although the bandwidth allocation methodology of the invention will be 
useful in any communications network connection between a sending and a 
20 receiving terminal, it is expected to find particular application in the allocation of 
bandwidth on Internet links with aggregate traffic, and such an application will be 
described hereafter as an illustrative embodiment of the invention. As a preface to 
the further description of the invention, the characteristics and parameters of such an 
Internet connection are first briefly considered. 
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The traffic on an Internet link usually arises from connections between pairs 
of computers. For each connection between a source and a destination computer, 
the information to be sent (typically a file comprising a number of bits or bytes) is 
broken up into packets at the source computer. The packets are sent from the source 
5 computer to the destination computer where they are reassembled to form the 
original file. 

In traversing the connection between the source and destination computers, 
the packets travel along a communications path across the Internet that consists of 
transmission links and routers. The source computer is connected to a first router by 

10 a transmission link; the first router is connected to a second router by another 

transmission link, and so forth. A router has input links and output links. When it 
receives a packet from one of its input links, it reads the destination address of the 
packet, determines which of the routers connected to it by output links should get 
the packet, and sends out the packet over the output link connected to that router. 

15 That process continues for each router in the connection with the final router 
sending the packet out on a link directed to the destination computer. 

At any given time, there are a number of simultaneous active connections 
using a given transmission link. The packets of the different connections are 
intermingled on the link; for example, if there are three active connections, the 

20 arrival order of 10 consecutive packets by connection number might be 1, 1, 2, 3, 1, 
1, 3, 3, 2, and 3. The packet intermingling is characterized as statistical 
multiplexing. Note also that links closer to either end of a connection are generally 
characterized by lower link speeds and fewer connections, while links traversing the 
network core are generally characterized by higher link speeds and more 

25 connections. 
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With the understanding that the number of connections using a given link 
varies over a time interval, a mean connection load over a stationary interval can be 
defined, which is designated herein as c. This mean connection load represents the 
average number of active connections using the link over a defined time interval. 

5 Another important traffic metric is the mean traffic bit rate over a stationary 

interval. That metric is determined by an evaluation of every packet arrival during a 
defined interval, including a determination of the size (in bits/bytes) for each 
arriving packet. The sum of the arriving packet sizes is then divided by the number 
of packets to arrive at a mean packet size (number of bits) in the interval. That 

10 mean packet size is multiplied by the packet rate for the interval to arrive at an 

average bit rate (bits/sec) on the link during that interval. That mean traffic bit rate 
is designated herein as x. (Alternatively, x can be computed as the sum of the 
packet sizes arriving in the interval divided by the interval length.) 

Note that as the number of active connections on a link increases, x tends to 

15 increase - i.e., as c increases, x tends to increase. 

When a packet arrives for transmission on a link, it is presented with a 
queuing mechanism, usually going into a buffer having a given capacity. From the 
buffer, the bits of the packet are put onto the link. The link has a fixed speed or 
bandwidth, p (bits/sec). The bits of a packet are put on the link at rate p. As will be 
20 apparent, the link bandwidth, p, must be larger than the mean traffic bit rate x. 
Accordingly, as x increases (often in association with an increase in the average 
number of connections on a link (c)) the link bandwidth will also generally need to 
increase. 
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If, when a packet arrives, other packets are already in the buffer awaiting 
transmission, then the arriving packet has to also wait for transmission. The time 
the packet must wait before it is actually transmitted on the link is characterized as 
packet delay and designated here as 8 (sec). Once the packet reaches the top of the 
5 queue (i.e., ready for transmission), the time required for transmission of the packet 
from the buffer onto the link is called the service time. 

Note also that, upon arrival of a new packet, if the buffer is full, the arriving 
packet is dropped. Both queuing delay and packet drops degrade the QoS of the 
traffic and one or more of these factors is often applied for defining acceptable QoS 
10 for a link. 

The QoS achieved on a link — particularly as related to delays and drops on 
the link — is closely linked to the traffic loading relative to the traffic capacity, or 
bandwidth, of the link. Stated differently, an increase in bandwidth for a given 
traffic load decreases service time, which reduces queuing delay and number of 
15 drops (and vice versa). Accordingly, poor QoS resulting from delays and drops on a 
link can be improved by increasing the link bandwidth. Note, as well, that loss and 
delay are also affected by the buffer size, and therefore that an increase in the buffer 
size will reduce the number of dropped packets. However, the queuing delay has the 
potential to increase with such an increase in buffer size. 

20 It should thus be apparent that a basic problem in bandwidth allocation is to 

choose the link bandwidth, /?, to satisfy the QoS criteria. The value of f5 so chosen is 
designated herein as the QoS bandwidth. A related parameter, utilization of a link 

(defined as: u = ) represents the fraction of the capacity of the link that is used 
on average. Since a router cannot transmit traffic faster than the link bandwidth, /?, 
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the link traffic rate, r, must be less than or equal to (3 in order to prevent major QoS 

problems. The QoS utilization then is the value of u determined as a function of the 
QoS bandwidth. 

The general bandwidth allocation problem can then be characterized as 

5 either: 

(1) What traffic load r can be put on a link with a specified bandwidth /? 
while maintaining a given QoS criteria; or 

(2) What link bandwidth /?is needed for a specified traffic load r in order to 
maintain the QoS criteria. 

10 For the illustrative embodiment of the invention described herein, the 

bandwidth allocation methodology of the invention is described in the context of 
Internet "best effort" traffic. It is noted, however, that the invention is intended to 
encompass an application of that methodology to other categories of traffic that are 
transmitted across a network. 

15 Note also that the heretofore-described delay and loss QoS criteria are 

applied in the development of the illustrative embodiment. It should, however, also 
be apparent that other QoS criteria could be used, and one skilled in the art will 
readily understand such alternate QoS applications from the description herein of 
the invention. 

20 For delay, two metrics are used: a target maximum delay, 5, and the delay 

probability, co (i.e., the probability that packet delay will exceed the target 
maximum delay). With respect to packet loss, it is assumed that a buffer size is 
chosen large enough that packet drops will be negligible - an assumption consistent 
with the current practice of service providers on their core links. To address the 
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problem of a large buffer size potentially leading to an increase in delay, appropriate 
limits are established for the maximum delay, <5, and delay probability, co in the 
formulation of the methodology of the invention. 

For purposes of the description of the invention herein, a packet stream 
5 should be understood to consist of a sequence of arriving packets, each with a size. 
The arrival number for the arriving packets will be denoted herein by v (v = 1 
signifying the first arriving packet, v = 2 the second packet, and so forth). The 
packet arrival times will then be designated by a v and the inter-arrival times 
designated by f v = a v+] - a v . Finally, the notation q y will be used to denote the size 
10 of the packet arriving at time a v . 

The QoS bandwidth for a packet stream depends on the queue-size time 
process, and the queue-size time process depends on the statistical properties of t v 
and q v . Accordingly, the statistical properties of the packet stream will be described 
herein by the statistical properties of t v and q v as time series in v. 

15 It is known that best-effort Internet traffic has long-range dependent t v and 

q v , but as the mean connection load, c, increases, the arrivals tend toward Poisson 
and the sizes tend toward independence. Relatedly, persistent, long-range 
dependent t v and q v have larger queue-size distributions than those for independent 

t v and q v (Citations to research in support of these propositions can be found in the 

20 companion cross-referenced paper for the reader wishing to pursue such). 

It is also known to those skilled in the art that long-range dependent traffic is 
burstier than independent traffic, resulting in smaller QoS utilization because more 
headroom is needed to allow for the bursts. It therefore follows that the bandwidth 
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allocation mechanism changes as the connection density, c, on a link changes 
because the statistical properties of t v and q v change - i.e., a "multiplexing gain" 
can be expected with higher utilization due to the reduction in dependence. 

The bandwidth allocation methodology of the invention encompasses a 
5 modeling approach base on queuing simulation, and the steps of that methodology 
are illustrated in the flow chart of Figure 1. As discussed above, the input to a 
queue is a packet stream defined by packet arrivals and sizes, and the methodology 
begins at step 101 with generation of a plurality of such streams (20 streams in the 
illustrative embodiment described herein). Each input stream has a traffic bit rate r, 

10 and a set of packet streams is chosen so that r varies over a wide range. As shown 
in step 102, many simulations are run for each stream, with systematically varying 
delays, S, and delay probabilities, co . For each simulation, the QoS bandwidth, (i , 
is determined for the fixed values of r, 5 and co . With the data from the 
simulations, and using known characteristics for queuing systems, a hierarchical 

15 statistical model is developed in step 103 (as described more particularly hereafter) 
that relates the QoS bandwidth to r, S and co . Finally, at step 104, the developed 
model is applied to determine bandwidth required to achieve a defined QoS criteria 
on a link. 

For a preferred embodiment of the invention, the packet streams are 
20 synthetically generated based on statistical models. Specifically, the synthetic 
packet streams are defined by arrivals and sizes generated by recently developed 
Fractional Sum Difference (FSD) statistical models for the t v and q v of best-effort 

traffic. A description of both the FSD model and the synthetic traffic streams 
developed therefrom is provided in Section 4 of the companion cross-referenced 
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paper, which also includes citations to additional research regarding such FSD 
models for the reader wishing to pursue additional detail. Inasmuch as the synthetic 
packet streams so generated are derived using known techniques and are used only 
as an input parameter for the modeling process of the method of the invention 
(described hereafter), detailed discussion of the FSD model is not warranted here. A 
brief overview of the discussion in Section 4 of the cross-referenced paper may, 
however, be useful. 

Synthetic packet streams generated by FSD statistical models have been 
found to provide an excellent fit for best-effort traffic, subject to the constraint that 
the number of connections per link (c) be above about 64 connections - a relatively 
small traffic load. 

Because the packet inter- arrival times, t v , and sizes, q v , are both non- 
gaussian and long-range dependent, a transform is made to a new time series in v 
that is close to gaussian. Thus t v and q v are each transformed to new time series, 

t v and q * , characterized as the sum of a gaussian long-range dependent time series 
and gaussian white noise. 

The marginal distribution of t v is modeled by a Weibull with shape X and 

scale a, a family with two unknown parameters. Estimates of X are almost always 
less than 1 . 

The marginal distribution of q v is modeled by a discrete-continuous 
distribution. The distribution is characterized by A atoms at sizes 0, (5) . . . (f> A s) , 
representing commonly occurring packet sizes. The corresponding atom 
probabilities are 0, (a) ... (f> A a) . The continuous part of the distribution is piecewise 
uniform on C intervals over the range of packet sizes. The breaks occur at C - 1 
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distinct values, $ b) . . . , greater than the packet size at and less than the 
packet size at (f>^ ) , and the endpoints of the intervals are these values together with 
the packet size at $ s) and the packet size at <fi ( A s) . The probabilities for the intervals 
are ^ (,) ... 0<°. 

A synthetic packet stream is generated by specifying the following 
parameters: (1) A atom probabilities ^ (a) at A sizes$ s) ,C-1 break points $ b) , and 
C interval probabilities $ l) ; (2) the packet rate per connection, y ; and (3) r . 

It is noted that, in the initial formulation of the methodology of the 
invention, the inventors also used live traffic data in their modeling process as a 
validity check on the use of the generated synthetic traffic streams. Each of the live 
streams of traffic constituted measurements of arrivals and sizes from Internet links. 
The data bases from which the live traffic was obtained and the relevant parameters 
for that traffic are described in Section 5.1 of the companion cross-referenced paper. 
From the inventors' work with the live and synthetic traffic streams, it became clear 
that the FSD synthetic streams provide a close approximation to measured statistical 
properties of the live data. There are also advantages to the use of the synthetic 
stream data over live traffic data: (1) the synthetic data allow extrapolation to traffic 
bit rates for which no live data are yet available; and (2) the synthetic data provide 
greater accuracy than the live traffic because one can generate arbitrarily large 
packet streams which may be needed to derive particular relationships (an 
impossibility with live data which are non-stationary through time and therefore 
must be broken up into individual packet streams over short periods of time). 

A synthetic trace database was constructed for application in the modeling 
methodology of the invention consisting of 20 traces, and each trace consisting of 

11 



Express Mail Label No. EU033990957US CAO 2-2-6 

20 million packets generated by the FSD statistical model. The 20 million packets- 
per-trace value was chosen so that the QoS utilizations for repeat runs of the 
simulation with the same synthetic traces and the same QoS criteria would differ by 
less than 0.005. 

As discussed above, the generation of a synthetic stream requires the 
specification of the marginal distribution of the packet sizes q v , the packet rate per 

connection y p , and the traffic bit rate r . In the illustrative embodiment described 

here those parameters are specified as follow. The size marginal distribution of all 
20 synthetic traces was taken to be the empirical marginal distribution of a set of 
live traffic traces which appeared to be representative. While one could also specify 
a model from the discrete-continuous class described above, this somewhat 
simplifying approach was follow because the result is very nearly the same. The 
packet rate per connection, y , was specified to be 1.5 packets/sec/connection; 

which was an average value experienced across all of the live traces evaluated by 
the inventors. 10 values of r were specified, with 2 traces at each of these values. 
The smallest value used for this illustrative embodiment was approximately 1 
megabit/sec, corresponding to the smallest bit rate in the live traces evaluated by the 
inventors. The 9 remaining values increase by successive powers of 2, so the largest 
value of r is about 512 megabits/sec. 

For each of the 20 synthetic traces generated as described above, each trace 
having a predetermined value of r , 25 runs were carried out for the illustrative 
embodiment, each with a particular queuing simulation. Specifically, for each run, a 
target delay £and a delay probability co were selected. Simulations were then run 
to find the QoS bandwidth fi — i.e., the bandwidth that results in delay probability 
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co for the delay 5. This also yields a QoS utilization u = rf (3 . In the illustrative 

embodiment 5 QoS delays (0.001 see, 0.005 sec, 0.010 sec, 50 ms, 0.100 sec) and 5 
QoS delay probabilities (0.001, 0.005, 0.01, 0.02, and 0.05) are used, and the 25 
combinations of the two QoS delay criteria were distributed among the 25 runs for a 
5 given trace. 

Each queuing simulation of a set of runs is characterized as follows: 

• the buffer is infinite; 

• the queuing discipline is first-in-first-out; 

• a trial QoS bandwidth is chosen; 

10 • the queuing delay at the arrival times of the packets is measured, which 

determines the queuing delay process; 

• from the queuing delay process the delay probability for the chosen S is 

found. 

The simulation is iteratively repeated, changing the trial bandwidth, until the 
15 attained delay probability approximately matches the chosen delay probability co . 
The optimization is straightforward because £ decreases as the trial bandwidth 
increases for fixed co . 

In the optimization, the utilization is capped at 0.97 - i.e., if the true QoS 
utilization is above 0.97, it is set to 0.97. The reason for this cap on the utilization 
20 value is that a logistic scale, log (u/(l - «)) , is used in the modeling of the data (as 
described below). Above about 0.97, the logistic scale becomes very sensitive to 
model mis-specification and the accuracy of the simulation, even though the 
utilizations above 0.97 are, for practical purposes, nearly equal. Similarly, the lower 
range of the utilization is limited to a minimum of 0.05. 



13 



Express Mail Label No. EU033990957US CAO 2-2-6 

The resulting data for each of the 20 synthetic streams includes 25 
measurements, one per run, of each of three queuing variables: the QoS utilization 
w, the delay S, and the delay probability co . For the described illustrative 
embodiment, the range of £is 0.001 sec to 0.100 sec, and the range of co is 0.001 to 
5 0.05. As already described, the range of the QoS utilization is 0.05 to 0.97. The 
additional variables that measure the statistical properties of the stream, such as r , 
are constant across the 25 runs for each stream. 

In the discussion hereafter of the data modeling according to the method of 
the invention, the following notation is used for the simulation data: S } for j = 1 to 5 

10 are the 5 values of the target delay and co k for k =1 to 5 are the 5 values of the delay 

probability. u ijk is the QoS utilization for delay Sj , delay probability co k , and 

stream i, where i = 1 to 20. 

An objective of the method of the invention is to find an algebraic 
relationship to address the dependence between QoS utilization, w, and the 
15 heretofore-discussed QoS criteria: delay 8 and the delay probability co . To that end, 
a general model for that dependency relationship was formulated of the form: 

where / is a function of u and g is a function of S and co that can vary from 
stream to stream, and the e iJk are independent random errors with mean 0. 

20 The inventors have determined that the function log 2 (S) provides a logical 

scale for the delay, S, and in an initial step for determining the dependency 
relationship, the conditional dependence of u on log 2 (S) , given co and the stream, 
is considered. 
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A powerful method for the consideration of such conditional dependence is a 
trellis display — i.e., for each combination of delay probability co k and trace i, plot 

u ijk against log 2 (Sj ) . Such a plot was developed for the 20 synthetic packet 

streams, which is shown in Figure 2 — utilization graphed against log delay 
( log 2 (S) ), given delay probability, co , and bit rate, r . As will be seen in the figure, 
there are 50 panels divided into 10 rows and 5 columns. On each panel, u ijk is 
graphed against log 2 (£ y ) for one value of co k and for two streams with the same 
value of T i . Note that the values of i and co k are shown in the strip labels at the top 
of the panel, and the value of r. for the two traces on the panel is given inside the 
panel. The t i are constant for all panels in a column, and increase as one moves left 

to right through the columns ~ in the leftmost column the value is about 1 
megabit/sec and increases by a factor of 2 in each successive column moving left to 
right. The co k are constant for all panels in a row, and increase as one moves from 

bottom to top through the rows. There are 10 points per panel, 5 for each of the two 
streams. Because the utilizations are so close for the two streams of each panel, each 
pair of points for a specific value of \og 2 (Sj) typically appears as a single point. 

A number of overall effects of r , S, and co on u can be seen from Figure 2. 
First, as noted above, there is little difference between the two sets of utilizations for 
each pair of streams with the same value of r ; in fact, the choice of a large sample 
size for the synthetic packet streams (20 million packets per stream), essentially 
ensured this result. Additionally, for each pair of values of co and r , there is an 
increase in u with increasing S, a strong main effect in the data. There is also an 
increase with r for fixed S and co , another strong main effect. A main effect for co 
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is also to be observed, but smaller in magnitude than for the other two variables. 
The dependence of u on log 2 (S) is nonlinear, and changes substantially with the 
value of r ; as r increases, the overall slope in u as a function of log 2 (£) first 
increases and then decreases. In other words, there is an interaction between 
5 log 2 (£) and r . 

From these effects observable in the trellis of Figure 2, the inventors have 
determined that a logit function would provide a good fit to the data for the function 
of u in the general model of the dependency relationship shown above. Thus 



(log 2 given delay probability, co , and bit rate,r , and using the same streams 

and layout as for Figure 2. As can be seen, the logit function does an excellent job 
of simplifying the dependence. The dependence on log 2 (5) is linear. There does not 
appear to be any remaining interaction among the three variables: log 2 (£) , r , and 

15 co . To help show this, 10 lines with different intercepts, but the same linear 

coefficient, have been drawn on the panels. The method of fitting will be described 
below. The lines provide an excellent fit. The only departure from the linear pattern 
occurs for the very highest values of r and the higher values of 8. However, the 
effect, while appearing large on the logistic scale is nearly inconsequential on the 

20 utilization scale since it involves values of u close to 97%, the maximum for that 
parameter. 

From consideration and extension of earlier work by Erlang in regard to 
delay probability (as described in Section 6.2 of the companion cross-referenced 
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Figure 3 provides a plot of logit utilization ( logit 2 (u ijk )) against log delay 
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paper), the inventors hypothesized that the function of delay probability, co , in the 
general model of the dependency relationship above, would be a term of the form 
-log 2 (-log 2 (co)) , the negative complementary log of co . The dependence of 
logit 2 (w) on -log 2 (-log 2 (a>)) was then evaluated by use of trellis plots similar in 
5 design to those of Figure 3. Figure 4 shows those plots for the 20 synthetic streams. 
On each panel, logit utilization (logit 2 (w l>Jt )) is graphed against the negative 

complementary log of the delay probability (-log 2 (-log 2 (a>)) ) for one value of Sj 
and for two streams with the same value of r i . As will be seen from Figure 4, the 
inventors' hypothesis is on target: logit 2 (w) is linear in -log 2 (-log 2 (&>)) , and the 

10 slope remains constant across streams and across different values of S. To help 
show this, lines with the same linear coefficient but different intercepts have been 
drawn on the panels. As with Figure 3, the method of fitting these lines will be 
described below. The lines provide an excellent fit except for the errant points for 
high utilizations observed earlier. 

15 The heretofore described model building and empirical findings (as shown in 

Figures 3 and 4) shows that a very simple model relates the QoS utilization u to the 
QoS delay, S, and the QoS delay probability, co . That model, which is hereafter 
characterized as the FSD delay model, is: 



20 where e is a random variable with mean 0 and variance cr 2 (e) that expresses 

random error variability of the model for predicting utilization. The parameter ju in 
the model changes with the packet stream being modeled and serves as a summary 




(1) 
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of the statistical properties of the stream. Hereafter, the dependence of /u on the 
measured statistical properties of a stream such as r will be considered. 

The model, rewritten as an equation that describes the specifics of the data 
here, is: 



log 



2 



u ijk 



= +o s log 2 {5 j )+o a (-log 2 (-log 2 (a k )))+€.., , (2) 



where the & ijk are independent realizations of the random variable e. 

Equation 2 was fitted to the data for the 20 synthetic packet streams in order 
to estimate the coefficients n,,o s and o a . Because of the outliers that occur at high 
utilizations, the fitting was carried out using the bisquare method of robust 
10 estimation. The estimates of o s and for the data are o 5 = 0.379 and o m = 0.863. 

The data fits are used in Figures 3 and 4, which show that the fitted equation lies 
very close to the data except for the small number of errant values. In other words, 
the FSD delay model provides an excellent fit to the data. 

To model the random variable e, the statistical properties of the residuals — 
15 the departures of the logit 2 (u iJk ) from the fitted equation - were considered. Except 

for a small fraction of errant residuals for high utilizations, e is well approximated 
by a normal distribution. The variance was estimated robustly by 0.75 times the 
median absolute residual - i.e., cr 2 (e) = 0.1 13. 

The coefficient \i in the FSD delay model of Equation 2 varies with the 
20 packet stream and reflects how the changing statistical properties of the streams 
affect the QoS utilization. Although the fact of a single numerical value |a 
characterizing how the statistical properties of a stream affect the QoS bandwidth is 
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part of the simplicity of the model, the variation with packet stream necessitates a 
secondary modeling of ja. Here, that modeling of ja as a function of packet stream 
statistical variables is considered. The end result is a hierarchical model for the QoS 
utilization: the first-level FSD delay model and a second-level model for \i. 

5 Based on deduction from the earlier referenced work of Erlang and 

consideration of changes in utilization, w, with change in connection density, r , as a 
function of both pure rate gain and multiplexing gain (as described in Section 7 of 
the companion cross-referenced paper), the inventors hypothesize that \x is linear in 
l°g2( r ) • F rom this hypothesis, the second level model for n is projected as: 

10 // = 0+o r (log 2 (r)-24)K (3) 

where ^is a random variable with mean 0 and variance <x 2 (0 that expresses 
random error variability of the model for predicting [i. The predictor k>g 2 (r) has 
been centered at 24 as the approximate center of the data values for the packet 
traces. The model, rewritten as an equation that describes the specifics of that data, 
15 is 

/ / / .=o+o r (log 2 (r | .)-24)+C (4) 
where the ^ are independent realizations of £ 

Figure 5 provides a plot of the estimated stream coefficients /u i of the fitted 
FSD delay model against log bit rates (log 2 (r g ) ) for the 20 synthetic packet streams. 
20 There are 20 points on the plot but the two results for each pair of streams with the 
same value of r f . are so close that they appear on the plot as a single point. The 
dependence of Jj ( on log 2 (r i ) is clearly linear, so the inventors' hypothesis is 
established. 
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The coefficients of Equation 3 were estimated by fitting it to the fi i and r i 

values using the bisquare method of robust estimation used above. The estimates for 
the synthetic packet streams are: 

o = 5.500 and 6 T =0.709. 

To model the random variable £ the statistical properties of the residuals — 
the departures of the fi i from the fitted equation, were evaluated. It was found that £ 

is well approximated by the normal distribution except a few errant traces with 
residuals large in absolute value. The standard deviation was estimated, as before, 
by 0.75 times the median absolute residual. That estimate is: 

&(£) = 0.036. 

The FSD delay model of Equation 2 (first-level model) together with the 
second-level model in Equation 4 for \x. (the stream coefficient model) together 
provide a hierarchical model for bandwidth allocation: predicting the QoS 
utilization, w, from 8, co and r . The model can be written in the form: 

logit 2 (u) = o+o r T+o s log 2 (<?)+<>„ (-log 2 (-log 2 , (5) 

where y/= e + C, is a normal random variable with mean 0 and variance 
a 2 (i//) = a 2 (e)+<r 2 that expresses random error variability of the model for 
predicting u. The model, rewritten as an equation that describes the specifics of the 
data here, is: 

togit^J^+o.r^ (6) 
The estimated values of the parameters from the synthetic packet trace data 

are: 
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6 = 5.500, 6 T = 0.709, 6 6 = 0.379 and 6„ = 0.863. 
The estimated value of cr{y/) from the data is: 
<x<y) =0.119. 

The hierarchical model here provides a model for predicting the QoS 
5 bandwidth given the delay, delay probability, and traffic bit rate. This is illustrated 
in Figure 6. On each panel, the predicted utilization, w, is graphed against log bit 
rate (log 2 (r)) from 1 megabit/sec to 1 gigabit/sec for the values of 8 and co shown 
in the strip labels at the top of the panel. There are six panels: two delay 
probabilities — 0.001 and 0.05 — and three delays — 0.001 sec, 0.010 sec, and 
10 0.100 sec. 

Herein, the inventors have disclosed a new and improved method for 
determining QoS utilization in an Internet link. Numerous modifications and 
alternative embodiments of the invention will be apparent to those skilled in the art 
in view of the foregoing description. In particular, it should be understood that, 
15 while the invention has been described in terms of best effort traffic over an Internet 
link, the principles of the invention will equally apply for other traffic types and 
other arrangements for transmitting data from a data source to a receiver. 

Accordingly, this description is to be construed as illustrative only and is for 
the purpose of teaching those skilled in the art the best mode of carrying out the 
20 invention and is not intended to illustrate all possible forms thereof. It is also 

understood that the words used are words of description, rather that limitation, and 
that details of the structure may be varied substantially without departing from the 
spirit of the invention and the exclusive use of all modifications which come within 
the scope of the appended claims is reserved. 
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